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BOUNDARY CONDITIONS ASSOCIATED WITH THE 
PAINLEVE III' AND V EVALUATIONS OF SOME RANDOM 
MATRIX AVERAGES 

P.J. FORRESTER AND N.S. WITTE 


Abstract. In a previous work a random matrix average for the Laguerre uni¬ 
tary ensemble, generalising the generating function for the probability that an 
interval (0, s) at the hard edge contains k eigenvalues, was evaluated in terms 
of a Painleve V transcendent in cr-form. However the boundary conditions for 
the corresponding differential equation were not specified for the full parameter 
space. Here this task is accomplished in general, and the obtained functional 
form is compared against the most general small s behaviour of the Painleve 
V equation in cr-form known from the work of Jimbo. An analogous study is 
carried out for the the hard edge scaling limit of the random matrix average, 
which we have previously evaluated in terms of a Painleve IIP transcendent in 
cr-form. An application of the latter result is given to the rapid evaluation of 
a Hankel determinant appearing in a recent work of Conrey, Rubinstein and 
Snaith relating to the derivative of the Riemann zeta function. 


1. Introduction 


The Laguerre unitary ensemble (LUEjv) refers to the eigenvalue probability 
density function (pxl.f.) 

1 N 

(i-i) ^n A “ e_Ai n 

N,a 1=1 1 <j<k<N 

where 

N-l 

(1-2) CW.a := I] r V + 1 ) r 0' + a + 1), 

j =o 


with support on A; € [0,oo). For a £ Z>o, this eigenvalue p.d.f. is realised by 
non-negative matrices X^X where X is an M x N complex Gaussian matrix and 
a = M — N. In our work [3] the average over (1.1) 


(1.3) E 2 ,N((0,s);a,fj,-,€) := 


C N , a 
C N , 


N 


a+fi 




l—l 


A 

/ LUEjv 
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where xlP = 1 f° r A i £ J, Xj ' 1 — 0 otherwise, and the normalisation is chosen so 
that 


(1.4) 


i?2,jv((0,s);a,/i;£) 


:= 1, 


s=0 


was characterised as a r-function for the Painleve V system. As a consequence, it 
was shown that 

(1.5) W N (s; a, /q £) := log ^"^^^((O, s); a, /q £)) , 
satisfies the Jimbo-Miwa-Okamoto cr-form of the Painleve V equation 

3 ■> 2 

(1.6) (scr'Jrf - 

with 


°V — SCTy + 2(o"y) 2 + (y ' 13j)t7y 
3=0 


+ 4 ff ('A + °v) — 0) 

l=o 


(1.7) z/o = 0, v\ = ~n, i >2 = N + a, 133 = N, ^ z/j = 2AT + a — q. 

j=o 

For this to uniquely characterise Wjy, a boundary condition must be specified. 
However in [3] only in the cases q = 0 and q = 2 were a boundary condition 
specified for general £. 

Also considered in [3] was the hard edge limiting average 

(1.8) i^ ard (f;a,/i;£) := lim ( CjV ’ a # 2 ,at((0, a, q; £)). 

N—>oo \C AT,a+Ai 4A / 


It was shown that 

(1.9) rphard 

where setting 

( 1 . 10 ) 


£r(t;o.wO=«P / u h {s-a^- 0 -, 

. n S 


u (s;a,q;£) = -(oiii'W 4-r-), 


the function <7m'(s) satisfies the Jimbo-Miwa-Okamoto cr-form of the Painleve 
IIP equation 

(1.11) (so-" n /) 2 — uir> 2 (o’n I /) 2 + cr( n /(4o-( n / — l)(oni' — soW) — ^(wi — V2) 2 = 0, 
with parameters 


(1.12) tq = a + q, V 2 = a — fx. 

Again only in the cases q = 0 and q = 2 were boundary conditions specified for 
general £. 

The aim of this work is to specify the boundary conditions relevant to both 
(1.5) and (1.6) for general values of the parameters. In the case of (1.5) this is done 
by writing the average (1.3) in its equivalent determinant form and evaluating the 
matrix elements in terms of certain \F\ hypergeometric functions. With the small 
s behaviour of the matrix elements determined, it turns out that the determinant is 
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such that its corresponding small s behaviour can readily be deduced. The small s 
asymptotic form of (1.3) then follows immediately. Scaling this asymptotic form as 
required by (1.8) then gives the small t behaviour of Ay ard (f: a, /x; £) and the small 
s behaviour of u h (s\ a, /z; £). 

The general small s asymptotic form of the permitted solutions of (1.6) and 
(1.11) have been given by Jirnbo [5]. As part of this study the boundary conditions 
found here are compared against these general forms. It is found that in both cases 
only one of the two branches permitted by the general solution is present in our 
random matrix problem. 

As an application of our results we specify the rapid computation of the power 
series expansion of a certain Hankel determinant of Bessel functions. The latter is 
known form our work [3] to be a special case of A’y ard (t; a. p; £). The coefficients 
in the power series appear in an asymptotic formula obtained recently by Conrey, 
Rubinstein and Snaith [1] for the integer moments of the derivative of the charac¬ 
teristic polynomial of a unitary random matrix. This in turn has application to the 
study of the derivative of the Riemann zeta function on the critical line. 

2. Small s Expansion of E 2 ,n{(0, s);a,/z;£) 

A standard result in random matrix theory, which in fact goes back to an identity 
of Heine (see [7]) expresses the random matrix average (1.3) as a determinant 

(2-1) E 2 ,n((0, s ); a, /x; £) = I ^' C ' N - a det [w j+k ]j,k=o,...,N-i, 

where 

pOO 

(2.2) w n := / dX (1 — £X(0, 5 ))(A - sf A"+“ e - A . 

Jo 

Unless fi is a non-negative integer (2.2) is not well defined for s real and positive, 
which is the domain of interest. To remedy this, we note that simple manipulation 
gives 

POO PS 

(2.3) w n := / dX (A - s)» X n+a e~ x + (1 - f) / dX (A - sfA n+a e -\ 

Js JO 

and in the second integral of this expression write (A — s) M = with 

—7T < arglog(A — s) < 7T. We then obtain 

POO PS 

(2.4) w n := / dX (A - sf A n+a e - A + (1 - fie™ 1 * / dX {s - X)> 1 X n+a e- x 1 

J s Jo 

which is well defined for 5ft(/i) > — 1 and 5ft(a) > — 1 for s > 0 with the additional 
constraint 5ft(/x + a) > — 1 at s = 0. 

We seek the leading terms in the small s expansion of (2.4). These can be 
read off from an explicit evaluation in terms of the \F\ confluent hypergeometric 
function [8]. 
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Proposition 2.1. Subject to the conditions $t(p) > —l, 3?(a) > —1, 3?(/x + a) > —1 
and p + a ^ Z>o we /lave 

(2.5) w n = a n (s) + s n+M+a+1 b n (s), 

where a n (s), b n (s) are analytic about s = 0 and given explicitly by 
a n {s) = T(p + n + a + \)e~ s \Fi{—a — n; —p — a — n; s), 


(2.6) b n {s) = 


T(p + l)r(n + a + 1) 

T(/r + n + a + 2) 


(1 — £)e’” M -: 


sm 7TO 


sin7r(/x + a) 


x e S i-Fi(/x+l;/x + fl + 7i + 2js). 


In particular, under the above conditions, 

(2.7) 


w n ~ a„(0) + sa(j(0) 

s—>0 

where 




a n (0) 

= T(p + n + a + 1), 

(2.8) 

<(0) 

= —pT(p + n + a), 


MO) 

r(/i + i)r(n + a + 1 ) / 

r(// + 71 + CL + 2) \ 


(i-£> 


sm 7ra 


sin7r(/x + a) 


Proof. The results (2.7) and (2.8) are immediate corollaries of (2.5) and (2.6) and 
the fact that 

1 ^ 1 ( 7 ; a; s) = 1 + — s + O(s 2 ). 
a 

To derive (2.5), we note that simple manipulation shows 

/•OO /»oo 

/ dX (A - s) A ‘A n+a e _A = s a+n e -5 / dA (1 + A/s) n+a A' i e" A . 

is io 


But with 


W k , m (z) = 


ykp-z/2 


T(l/2 - k + m) J o 
specifying the Whittaker function, it is known that [8] 

T(2m) 


/*oo 

/ dt (1 + t/z) k -^ 2 +m t- k - l/2+"* e -*, 

Jo 

Mk-m(z) 


where 


T(—2m) 

= T(l/2 — fc — m) Mk,m( z ) + r(1/2 _ k + m) 
Mk, m (z) = 2 m+1/2 e" z/2 iFi(l /2 - fc + m; 2m + 1; 2). 


Consequently 

/•OO 

(2.9) / dA (A — s) A1 A" +a e -A = T(/i + a + n + l)e _s iJ 7 ’i(—a — n; — p — a — n; s) 

J S 

+ r(M + ~ Q ~ n ~ + 1; /x + a + n + 2; s). 

1(—a — n) 

The left-hand side of (2.9) exists for 3J(/z) > — 1 if s > 0 and 3 ?(p + a) > — 1 if 
s = 0, whereas the right-hand side is valid in this parameter domain except for 
p + a + n £ Z> 0 , and in this case the individual terms have a simple pole at 
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a + n (j Z>o or are undefined when a + n £ Z>o- Needless to say the sum of the 
terms on the right-hand side has the same analytic character as the left-hand side. 

Regarding the second integral in (2.4), we first note that a simple change of 
variables gives 

f dX (s — \y\ n+a e- x = S n + l + a +^ e ~ s f dx (1 - x) n+a x^e sx . 

Jo Jo 

But the integral on the right hand side is the Euler integral representation of the 
i Fj function, which shows 


( 2 . 10 ) 



\yx n+a e~ x 


r(/x + i)r(a + n +1) 

r(/x + a + n + 2) 


s v. +a +n+i e -s iFi ^ + l-n + a + n + 2 -s). 


This latter relation is valid for 3i(/n) > —1 and 5ft(a) > —1 when s > 0. Substituting 
(2.9) and (2.10) in (2.4) and using the appropriate gamma function identities gives 
(2.5), (2.6). □ 


When /r + a £ Z>o we have to consider two exceptional cases where one of the 
hypergeometric functions are not defined - the first when a + n € Z>o for which the 
hypergeometric function is indeterminate, and the second when a + n ^ Z>q and 
the hypergeometric function has a simple pole. These two cases can be recovered 
by taking suitable limits and we just state the final results. 

Proposition 2.2. When fj, + a = j £ Z>o and a + n = k £ Z>o with n + j > k we 
have 

+ (- l) n+J+k (1-0 i n T J Tnl sU+j+1 r^i (n + j + 1 - fc; n + j + 2 ; s) l , 

(n + j + l)! J 

and to leading order in small s we have 

( 2 . 12 ) 

w n ~ (n + j)\ -(n + j-k)(n+j-l)\a + (-l)^ +fc (l - 0 ^±lZ^L s n+5+\ 
s^o (n + j + 1)! 

Note that the condition n + j > k is the same as fi > 0, which falls within the 
domain of interest. The key difference of (2.12) with (2.7) and (2.8) is that the 
non-analytic term is now polynomial and the second part of this term is absent 
having been cancelled by a counterbalancing term. 
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Proposition 2.3. When fi + a = j £ Z>o and a + n (f Z>o we have 


(2.13) 


' n+j 

E 


(-a - n)i(n + j - l)\ t 


+ 


K 1=0 

r(ti + l)T(a + n + 1) 
(n + j + 1)! 


l\ 


(1 - e)e iwM s n+J+1 iPi(M + 1; n + j + 2; a) 


, , lU - sin7rar(^ +l)r(g + rt + 1) n+j+1 
[ ’ 7r (n + j + 1)! 

x E + 1 )+i > (n + j + l + 2) - tjj(n + l + 1) - logs] yrl, 

^ (n + j + 2), l\ j 


and its leading order behaviour for small s is 


(2.14) ~ (n + j)! + (a-j)(n+j-l)!s 

s —^0 


. ( a j)w+j + l n +j+l 

(n + j + 1)! 


7re Z7ro 
sin 7ra 


(1 — 0 + V'(l) + V*(n + j + 2) - ^(/x + 1) - log s 


The expansion (2.14) differs significantly from (2.7) and (2.8) because of the 
presence of logarithmic terms which now replace the non-analytic contributions of 
the generic case. 


Corollary 2.1. Under generic conditions on /i + a we have 


(2.15) det['U)j + fe] Ji fe = o,...,iv-i 

= det[T(/x + a + 1 + j + k)]j t k= o,...,iv-i 
— //sdet[r(/i + a + j) T(/i + a + 1 + j + fc)] 3 =o,...,n-i + 0 (s 2 ) 

fe=l,... ,N — 1 

+ s 7‘+“+ 1 5 0 (o) det[T(p, + a + 3 + j + &)]j,k=o,...,iv —2 {1 + O(s)} 

+ 0(s 2(/i+a+1) ). 

Proof. According to (2.7) 


det [wj + k]j,k=o,...,N-i 

a ~ o det[a i+fc (0) + sa' j+k ( 0) + s^ +a+1+i+k b j+k { 0)\ jlk =o,...,N-i 
~ det[oj + fc(0)]j |fc= o,...,iv-i + 5[s] det [iij-i-fc(0) + 5 (Xj_j_^.(0)]j ) / c= o,... ) iv— l 

+ s /i+a+1 6 o (0) det[a i+fe+2 (0)] ijfe= o,...,Ar-2, 

where [s]P(s) denotes the coefficient of s in P(s). Recalling the explicit formula 
for a ra (0) as given in (2.8) we obtain the constant term and the term proportional 
to gA‘+ a + 1 in (2.15). It remains to compute the coefficient of s, which according to 
(2.8) has the explicit form 


(2.16) 


[s] det[T(/i + a + 1 + j + k) - jj,sT{n + a + j + fc)kfc=o,...,iv-i- 
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Using the linearity formula 

det[ai • • • aj + bj • • • a„] = det[ai • • • a j ■ ■ ■ a„] + det[ai • • • b j ■ ■ ■ a„], 

where the a’s and b’s are column vectors, on each column of the determinant we 
see that of the terms proportional to s only the one obtained from expanding the 
first column in non-zero (all the rest result in two identical columns), and the 
determinant given by (2.15) results. □ 


It remains to evaluate the determinants. For this task we make use of the 
identity [6] 

n —1 

det [r(zfc + j)]j,fe=o,...,n-i = r(^) (z k - Zj). 

j —0 0<j<k<n —1 

After straightforward manipulations, gamma function evaluations of all the deter¬ 
minants in (2.15) can be obtained. Substituting in (2.1), and recalling that the 
normalisation is such that at s = 0 En is equal to unity, we obtain the sought 
small s expansion of Ejy and thus Wn valid for general values of the parameters. 


Proposition 2.4. For IR(yi) > —1, 3J(a) > —1 and fj, + a ^ Z>o we have 


(2.17) E 2 ,n(( 0, s); a, /x; 0 = 1 -:— s + 0(s 2 ) 

H + a 

r(/i + i)r(a + i)r(/x + a + n + i) /, ,, 
r 2 (A/ + a + 2)r(Ai + a+l)r(A0 V * 


sin 7T a 


siii7r(Ai + a) 


s^ +a+ i {1 + O(s)} 


+ 0 (. 


.2(li+ a +I) 


), 


and consequently 


(2.18) Wjv(s;a, A 4 ;?) = —Nfj, -s + 0(s 2 ) 

A i + a 

r(/i + l)r(a + + a + N + 1) ( _ . ifl _ sin 7ra \ 

r(Ai + a + 2 )r 2 (At + a + 1)T(N) \ siii 7 r(At + a) J 


s M+ a +i ^ + 0 (s)} 


+ 0(s 2( ' i+a+1) ). 


In the first exceptional case fi + a = j G Z>o and a = k G Z>o with j > k 
one can still use (2.17) but omitting the term involving the ratio of sines, in the 
case j = 0, or the whole term if j > 0. The situation of the other exceptional case 
A i + a = j G Z>o and a ^ Z>o is more complicated and more so for larger j , and we 
only give the examples of j = 0,1. 


Proposition 2.5. For > —1, 5ft(a) > —1 with n + a = 0 we have 


(2.19) E 2 ,jv((0, s);a,At = -a;£) = 1 

+ { - 1+ Ka e~” 1 ’ a (l-£) + a[2?/>(2) + ^(1) - tp(l - a) - ip(N + 1) - logs] \Ns 
l sin ira ) 

+ 0 (s). 
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For fx + a = 1 we have 

(2.20) ^2,n((0, s); a, /x = 1 — a; £) = 1 + (a — 1 )Ns 

+ a ^~ ^ {-r^e- i7va (l-0+2Tp(3)+'il>(2)-Tp(2-a)-i>{N+2)-logs}(N+l)Ns 2 

+ o (s 2 ). 


3. Comparison with the Jimbo solution 

The small s expansion of the most general solution permitted by (1.6), or more 
precisely its corresponding r-function (see (3.2) below) has been determined by 
Jimbo [5]. However in [5] the equation (1.6) is not treated directly. Instead the 
discussion is based on the equation 


(3.1) «") 2 - [C - s? + 2(C) 2 - (20 o + MC'] 2 

+ 4C , (C / — 0o)(C' ~ 2 — 0 S + 0oo))(C — — (0o + 0« + 0oo)) = 0, 

and the small s behaviour of the corresponding r-function ry(s), specified by the 
the requirement that 

(3.2) C( s ) = s ~r log Ty(s) + -(0o + 0oo)s + -j- [(0o + 0oo) 2 — 0 2 ], 

as 2 4 

was determined. 

Comparison of (3.1), (3.2) with (1.6), (1.5) shows that for the parameters (1.7) 

(3.3) £ 2 , iv((0, s); a, W £) = S ^+N(a+^) e -{N+a/2)s Tv ( s)) 
while in general 

(3.4) 0 O = —Zd, 0 s = i / 2-i / 3, 0oo = vi - v 2 - v z . 

Note that for the parameters (1.7) we thus thus have 


(3.5) 


#0 1-^1 @s — @00 


2N — a — [i. 




BOUNDARY CONDITIONS ASSOCIATED WITH THE PAINLEVE III' AND V... 


9 


The relevant result from [5] can now be presented. It states that the most 
general small s behaviour of ry(s) permitted by the equation (3.1) is 


(3.6) r v (s) = Cs<S -O' 4 /1 - M<>; ~ f + <t2) . 


r 2(_ CT ) r(i + ds + + g )r(i + 9s ~ ? 2 ° + g )r(i + e °° + a ) 

r 2 (2 + cr) r ( 9 s + 9 o~a2° ~ a )r( ^°°,^ ° ) 

9 S + 00 - CbTV-l I fls — 0Q — , $no ~ OA 


,1 + CT 


+ 


i r 2 (a) r(i + t 7 ^^- g )r(i + 


-)r(i 


u r 2 (2 — cr) 


r( 


+ ffn + £ 


)T(f 




+ 0(|s| 2(1 - K(ff)) ) 


where C is a normalisation constant, while u and cr are arbitrary parameters. The 
above result was derived subject to the conditions 9q : 9 s ^ Z, i( 9oo ± cr) ^ Z, 
5 (0 s ±0o±cr) ^ Z and that 0 < 3J(tr) < 1 (a distinct solution was presented for a = 
0). These conditions therefore strictly apply only to the generic or transcendental 
solutions of the fifth Painleve equation. For generic parameter values the terms 
given in (3.6) uniquely specify all the subsequent terms in the convergent Puisuex- 
type expansion for £(s) about s = 0 

OO 

(3.7) c w = EE c ^ +b 

J=0 \k\<j 

i.e. with any two of cycpCpi or cy-i given. 

To relate this to E 2 ,n, we see from (3.3) and (3.5) that we require cr 2 = (a + n) 2 
and thus we can choose 


(3.8) 


cr = a + fJ,. 


This relation, cr = 0 O + 0 S , is a violation of one of the strict conditions given above 
and is in fact a sufficient condition for a classical solution, along with the necessary 
condition 9q + 9 S + 9^ = — 21V £ Z, which is the type of solution that we are 
dealing with here. However we conjecture that Jimbo’s conditions can be relaxed 
to accommodate such solutions and the corresponding formulae (or limiting forms 
if necessary) still hold. With this choice of a the coefficient of s 1 ~ a in (3.6) contains 
a factor of 


1 


r( 


1 

rpv) 


7 
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and thus vanishes. Simplifying the other terms gives 


t v (s) ~ Cs- N2 ~ N ^ { 1 + (2i !, + Q + M) % 


2 (a + /i) 

sin7r/i r(a + l)r(/x + l)r(iV + 1 + a + /x) 1+a+ „l 
U S in 7r(a + p) r 2 (2 + a + p)T{l + a + /x)r(AV) S ]' 


Substituting in (3.3) we see that this is in precise agreement with (2.17) provided 
we choose 


(3.9) 


u 


sin 7r/i 
sin7r(a + /x) 


(1 - £)e™ M 


sin 7ra 

sin7r(a + /x) 


4. The hard edge limit 

The hard edge limit is defined by (1.8). However, only in the cases p = 0, p = 2 
do we know how to prove its existence for general £ (in the case /x = 0 E^n can 
be written as a Fredholm determinant, while the case p = 2 is related to this via 
differentiation). However a log-gas viewpoint ([2]) indicates that the limit will be 
well defined, and moreover we expect that it can be taken term-by-term in the small 
s expansion of E^.n ■ In this section we will show that taking the hard edge limit 
of the small s expansion (2.17) give rise to an initial condition for the solution of 
(1.11) consistent with that allowed by Jimbo’s theory of the small s expansion of the 
Painleve III' equation. From a practical perspective this specifies i? 2 ard for general 
values of the parameters according to (1.9), while from a theoretical viewpoint it 
lends weight to the belief that (1.9) is indeed the correct limiting evaluation for 
general values of the parameters. 

Under the assumption that the hard edge limit can be taken term-by-term in 
the small s expansion of Proposition (2.4) is immediate. 


Corollary 4.1. For 3J(/x) > — 1, 9?(a) > — 1 and p + a ^ Z>o we have 


(4.1) E^(s;a,p-,0 = 1- 

r(/x + l)r(a+l) 
r 2 (/x + a + 2 )T(p + a + 1) 


A* 


4(a + /x) 


s + 0(s 2 ) 


(1 - t)e^ - - 


sm 7ra 


sin 7 r(/x + a) 


/ S \ M _ r a_ r 1 

( 4 ) f i+0 <»» 

+ 0(s 2(/i+a+1) ), 


and consequently the cr-function crin'(s) in (1-10) has the small s expansion 
ix(n + a) p 


(4.2) (Tin' (a) - - + 4( ^ a) a + 0(s 2 ) 

r(/x + l)r(a + 1) ^ sm7ra 


r(/x + a + 2)r 2 (/x + a + 1) 


(1 - £)e mtl - 


sin 7 t(h + a) 


/ s \ 

(i) {‘ + 0(<» 

+ 0(s 2(/i+a+1) ). 
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Some examples of exceptional cases not covered by the preceding corollary are 
the following. They are obtained by taking the hard edge limit of (2.19) and (2.20). 

Corollary 4.2. For 3?(/x) > — 1, 3?(a) > —1 and n + a = 0 we have 
(4.3) E 2 hard (s;a,/x = -a;e) = l 

+ { - 1 + -r——e~ 7ria (l -0 + a[2ip{2) + V>(1) - tHl - a) - log(s/4)]U 
l sin7ra J 4 

+ o(s), 


whilst for fj, + a = 1 we have 


(4.4) E% aid (s;a,iJ, = 1 - a; £) = 1 + (a - 1)| 

+ -0+ 2^(3) + m - W - a) - log(s/4)} (|) 2 

+ o(s 2 ). 

To compare these results to the small independent variable expansions given by 
Jimbo in the theory of III' , we must first undertake some preliminary calculations 
as the equation (1.11) is not directly studied in [5]. Rather the equation studied is 

(4.5) (to 2 = 4C'(C' - 1)(C - to + - nc)\ 


where we have identified do = — V 2 , #oo = ~v\ (9 q, Ooo are the parameters appearing 
in [5]). In terms of f(t) the r-function rnr(t) is specified by the requirement that 

(4.6) ((t)=t^ log mi' (t) + V2 4 Vl + t, 

and it is the small t expansion of Tm'(f) presented in [5]. Comparison of (4.5) and 
(1.11) shows that 

., „. . , , vAvo — Vi) s s 

(4.7) CM = -o'hf(s) H -- b i= 4' 

Recalling (1.10), (1.9), (4.7) and (4.6) we see 


( 4 . 8 ) 


£“( S ;a, /1 ;()=^M/ 4 T IIP ((). 
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In [5] the most general small t expansion of tiii/(£) as permitted by (4.5) is 
presented. It reads 


(4.9) mv(t) = c& 2 ~^/ A ji + 

r 2 (—g) r( 1 + a+^)r( 1 + ^) 

U r 2 (2 + a) T (V2 ^JL) T (UJl^JL) 

1 r 2 (cr) r(i + ^iZ)r(i + 
u r 2 (2 — a) r(^+^)r(^^) 


+ 0(|t| 2(1_SR(ff)) ) 


where as in (3.6) C is a normalisation, while u and er are arbitrary parameters. This 
result was established under the assumptions that 1 (y\ ± a) ^ Z and i (i >2 ± <r) ^ Z 
along with 0 < Sifter < 1 (for a = 0 a distinct solution is given). 

To see that this structure is consistent with (4.1) and (4.8), recalling (1.12) we 
see that for the right hand side of (4.8) to tend to 1 as t tends to zero we must have 
C = 1 and er = ±ui. Again this is a violation of first condition given above but we 
conjecture that the formulae have meaning under the following limiting procedure 
and are correct. Choosing the positive sign for definiteness, and then writing 

u u(v i — a) 

T( vi^_g_) 2T(1 + ^^) 

we see that requiring 

u. . sin nvi 

— (ui — a) —> u - as a— 

2 7r 

(4.9) reads 


(4.10) + 

_sin7r(ui - u 2 )/2 T(1 - l ' 2 % ~ 1 )r(l + V2 ~^ Vl ) / A 1+1,1 
U sin7rui T 2 (2 + ui)T(l + vi) \4y 

Recalling again (1.12) and (4.8) we see that this agrees with (4.1) provided 


(4.11) 


Sill ITfJ, 


siii7r(a + n) 


= (i-0< 


sm 7ra 


sin7r(a + p.) ’ 


(cf. (3.9)). 
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5. Application 

In a recent work relating to the application of random matrix theory to the study 
of moments of the derivative of the Riemann zeta-function, Conrey, Rubinstein and 
Snaith [1] obtained two asymptotic expressions associated with the derivative of 
characteristic polynomials for random unitary matrices. With U a Haar distributed 
element of the unitary group U(N), and e 101 ,..., e l0N its eigenvalues, let 

N 

(5.1) A A ( S ) = I]( 1 -se-^), 

3=1 

and 

(5.2) Z A {s) = e-^/V^i e " /2 s - JV/2 A A (s), 

(note that Z A (e“) is real for 6 real). In terms of this notation, the two results from 
[1] are 

(5-3) (|A^(l)| 2fc ) Ae u(iv) ~ b k N k2 + 2k , 

—*oo 

where 

(5.4) 6 fc = (-l) fc ( fc+1 )/ 2 Wj"V + /i)! 

h =0 ' ' 

x [x k+h ] [e~ x x~ kl ^ det[/ a+ /3_i(2v / )r)]a,/3=i,...,fc^ , 

and 

(5-5) (\Z' A (l)\ 2k ) A eu { N) ~ b’ k N k2+2k , 

rJ —KX> 

where 

(5.6) b' k = (_i)fc(fc+ 1 )/ 2 ( 2 k )\[x 2k ] ( e-*l 2 x- k2 l 2 det[/ a+/3 _ 1 (2^)]a, /3 =L..., fc ) • 

In (5.4) and (5.6) the notation [x p ]f(x) denotes the coefficient of x p in f(x). 

The relevance of these formulae to the present study is that the determinant 
therein can be identified in terms of E 2 ard . Thus, we have shown in a previous 
study [3] that for a € Z> 0 

( n 

J e~ s/A det[/ /i+a _ /3 (v / s)]a > /3=i,...,a. 

where 

m \ i TT 0 + M — I)- 
A(a,/r) = a! j j- - -. 

3=1 J ' 


(5.8) 
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Interchanging row /? by row a — (3 + 1 (/? = 1,..., a in order) we see from this that 


bk = 


k 


(5.9) 


A(M)^ o 


(k + h)\[x k+h ]EF xA (4x\k,k\Z = 1 ) 


v - C" 1 )* 

fc A(fc,fc) 


(2fc)![x 2fc ] (e x/2 ^ ard (4x; fc, fc; £ = 1)) 


Note that the Painleve III' parameters appearing in this solution are p, = a = 
fc € N and n + a = 2k £ 2N and thus we are dealing with the exceptional case 
of indeterminacy referred to in Section 2. However as was noted there the generic 
formulae still hold with to the modifications discussed and in particular the a- 
function has a small argument expansion of a purely analytic form. 

From ([4]) it is known that the determinants in (5.4) and (5.6) can also be 
expressed as a particular generalised hypergeometric function. Such an observation 
implies, for instance, that 
(5.10) 


x fc2 / 2 det[/ a+j g_i(2 v / *)]a,/3=i,...,fc = I 


3'- 


=i r (j + fc) 


F^ 1] (;2k-,xi,... ,x k )\ Xj =x 


where c; x \,..., Xk) has a series development about X\.,Xk =0 with an 

explicitly given coefficient for an arbitrary term. However this is not a practical 
or efficient way to compute the coefficients required in (5.4) or (5.6) for moderate 
or large fc as it involves the hook lengths of Young diagrams associated with the 
partitions of fc. 

According to (1.9), (1.11) and (4.2) 

( pAx 7 

-J — (oiif(s) + fc 2 ) 

where <Jm> (s) satisfies the particular er-Painleve IIP equation 

(5.12) ( so iii') 2 + (Tjjjf (4<r( n / — 1 )(ctih' — ser( n /) — — = 0, 

subject to the boundary condition 

(5.13) criir(s) ~ -fc 2 + ^ + 0(s 2 ), fc £ N. 

s—>o 8 

Substituting 

(5.14) o-iif(s) =f?(s) + 

(5.12) reads 

(5.15) (sr/") 2 + 4((r/') 2 - ^-)(r/ - sr/') - ^ = 0. 

We see immediately that r/(s) can be expanded in an even function of s about s = 0, 


7 ?( s ) = ^2 c nS 2n , 

n =0 


(5.16) 


cq = — fc 2 , fc £ N. 
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Moreover the coefficients can be computed by a recurrence relation. 
Proposition 5.1. Substituting (5.16) in (5.15) shows 


(5.17) 


1 

64(4fc 2 - 1) ’ 


while for p > 2 


(5.18) 


1 


2c\p{2p — 1) + (2 p — 1)/64 — 8pk 2 Ci 
/ p-2 

x 4fc2 + 1 )(p _ l)ci+lC p _i 


1=1 


p-2 


— ^ 2(1 + l)(p — l)(2l + 1)(2 P — 21 + l)c; + iCp_i 


1=1 


P-1 


- ^(1 — 2l)ciA p _i_ 1 ), 


i=i 

where 

q 

(5.19) A q = ^(l + 1 )(q - l + l)Cl +1 C q -l +1 . 

1=0 

Proof. With hi := (/ + 1)(2Z + l)c; + i we see 

oo p 

(5.20) (st/') 2 = 4 J2 h p-is 2p , H p = J2hih P -i 


p=i 


1=0 


and similarly with ai := (/ + l)c/_|_i we have 


(s?/) 2 = 4s 2 ^2 A p s2P i A p = y ^aiap-i. 
p= o ;=o 


It follows from this latter result that 
(5.21) 


(V) 2 “ “ S??,) = £ Gps2P ’ 

' ' p =o 


where 

p 1 

G p - 2l)cib p -u b 0 = - —, b p =4A p _ 1 (p> 1). 

1=0 

Substituting (5.20) and (5.21) in (5.15) and equating like coefficients of s 2p to zero 
shows that for p > 1 

H p — i + G p = 0. 

This for p = 1 implies (5.17), and for p > 1 implies (5.18). □ 
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Using Proposition 5.1 it is straightforward to calculate, via computer algebra, 
the first k coefficients in (5.16) for any particular value of k. Furthermore use of 
computer algebra gives the power series up to x 2k of 

^ ard (4x; k, k; £ = 1) and e x/2 E% ard {4x; k, k; £ = 1), 


according to (5.11). From these power series the formulae (5.9) are used to compute 
bk and b' k . In [1] the first 15 values of both bk and b' k were tabulated. This can be 
rapidly extended using the present method. However the resulting rational numbers 
quickly become unwieldy to record. Let us then be content by presenting just the 
16th member of the sequences, 


bi6 = 


^16 — 


307-23581-92867-760550281759 

2 272 gl30. 5 66. 7 42 n 24 13 21 17 ie. lg 14 23 10. 2 g6 31 5 37 3. 41 2 4 g2 47 gg 5g gl ’ 

4148297603-7623077808870586151748455369217213506671334530597 

2 264. 3 133. 5 66. 7 42 1:L 25. 13 21 17 16. lg 14. 23 ll. 2g 7. 31 6. 37 3. 41 2. 43 2. 47 . 53 5g . 61 


6. Acknowledgements 

This work was supported by the Australian Research Council. PJF thanks M. 
Rubinstein for relating the results of [1] before publication, and for the organisers 
of the Newton Institute program ‘Random matrix approaches in number theory’ 
held in the first half of 2004 for making this possible. NSW wishes to thank the 
organisers of the CRM program ‘Random Matrices, Random Processes and Inte- 
grable Systems’ held in Montreal 2005 for the opportunity to attend the program 
and in particular the hospitality of John Harnad during his stay. 

References 

[1] J.B. Conrey, M.O. Rubinstein, and N.C. Snaith. Moments of the derivative of the riemann 
zeta-function and of characteristic polynomials, math.NT/0508378, 2005. 

[2] P. J. Forrester. Log Gases and Random Matrices. 
http: //www.ms .unimelb. edu. au/~matpjf /matpjf .html. 

[3] P. J. Forrester and N. S. Witte. Application of the r-function theory of Painleve equations to 
random matrices: PV, Pill, the LUE, JUE and CUE. Commun. Pure Appl. Math., 55:679-727, 
2002 . 

[4] P. J. Forrester and N. S. Witte. Discrete Painleve Equations, Orthogonal Polynomials on the 
Unit Circle and n-recurrences for averages over U(N) - Pni' and Pv T-functions. Int. Math. 
Res. Not., 2004(4): 159-183, 2004. 

[5] M. Jimbo. Monodromy problem and the boundary condition for some Painleve equations. 
Publ. Res. Inst. Math. Sci., 18(3):1137-1161, 1982. 

[6] J.-M. Normand. Calculation of some determinants using the s-shifted factorial. J. Phys. A, 
37(22):5737-5762, 2004. 

[7] G. Szego. Orthogonal Polynomials. Colloquium Publications 23. American Mathematical So¬ 
ciety, Providence, Rhode Island, third edition, 1967. 

[8] E. T. Whittaker and G. N. Watson. A Course of Modern Analysis. Cambridge University 
Press, Cambridge, 2nd edition, 1965. 




BOUNDARY CONDITIONS ASSOCIATED WITH THE PAINLEVE III' AND V... 


17 


Department of Mathematics and Statistics, University of Melbourne,Victoria 3010, 
Australia 

E-mail address: p.forrester@ms.unimelb.edu.au 
E-mail address: n.witte@ms.unimelb.edu.au 


